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1. Introduction 

One of the major goals of lattice QCD is a precise determination of hadron form factors. 
Electromagnetic (EM) and scalar form factors of pions and kaons are fundamental observables in 
hadron physics. A detailed comparison between their chiral behavior on the lattice and that in 
chiral perturbation theory (ChPT) may provide a good testing ground for these theoretical tools as 
well as a determination of unknown parameters in ChPT, namely the low-energy constants (LECs). 
Reliable calculation of the form factors of semileptonic weak decays, such as the K^nlv decays, 
is important in the search for new physics through a precise determination of CKM matrix elements. 

In this article, we report on our calculation of these light meson form factors. Chiral symmetry 
is exactly preserved in our simulations, that enables us to directly compare our results with ChPT. 
We use the all-to-all quark propagator in order to accurately calculate connected and disconnected 
correlators of pions and kaons. 

2. Numerical simulations 

Our gauge ensembles of Nf=2+\ QCD are generated using the Iwasaki gauge and the overlap 
quark actions with a topology fixing term |[l]], which remarkably reduces the computational cost. 
While we only explore the trivial topological sector at the moment, the effect of fixing topology on 
the pion form factors turned out to be under control, typically a few %, in our previous study in two- 
flavor QCD [§]. On a A^.^' xiV, = 16^ x 48 lattice, we simulate a single lattice spacing a=0. 112(1) fm 
leaving a quantitative estimate of discretization errors for future studies. Four values of degenerate 
up and down quark masses m„d = 0.015, 0.025, 0.035 and 0.050 are taken to explore a range of 
the pion mass 310<M;j[GeV] < 560. Measurements of the form factors are carried out with the 
periodic boundary condition at a single value of the strange quark mass = 0.080, which is very 
close to the physical value 0.081 fixed from Mk- We have accumulated 2,500 HMC trajectories at 
each combination of m„j and m^. 

We calculate meson correlators using the all-to-all quark propagator [||]. Let us consider an 
expansion of the propagator using the eigenmodes {Xk-, u^) {k = \,.. . , llN^Nt) of the overlap-Dirac 
operator, namely D^^{x,y) = T.k^k^Uk{^)'^l(y)- It is expected that low-lying modes dominantly 
contribute to low-energy observables, such as the form factors. We exactly evaluate this important 
contribution by using 160 eigenmodes for each configuration. Remaining contribution from the 
higher modes is estimated stochastically by using the noise method. We refer readers to Ref. [Q] 
for further details on our measurement method using the all-to-all propagator. 

We calculate three- and two-point functions 

^^tiys x,tx',x" 

^^-!P'(X"-X')-!P(X'-X)^ (2.1) 

c'4:{At,p) = _L££(p;,(x',f+AOp;(x,o>e'-*('''-''), (2.2) 

where (P'^\) represents an interpolating operator for the initial (final) meson with a smearing 
function <p {(p'), and is either the EM current (7^), weak current (V^), or scalar operator (S). Using 
the all-to-all propagator, we can accurately calculate these correlators by taking the average over 
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Figure 1: Left panel: effective value F^^(Af , Af';^^) with local (open symbols) and smeared operators for 
pion source and sink (filled symbols). Note that we can take arbitrary combinations of (Af , Af') by using the 
all-to-all propagator. Right panel : Fy^{q^) as a function of q^. We also plot the q^ dependence expected 
from the VMD hypothesis by the dashed line. 



the location of the meson source, namely (x,t) in Eqs. (2J^) and (22). In addition, the correlators 
with different choices of ^^'^ and p^'' can be calculated with small additional costs. In this study, we 
take two choices of ^^'\ namely local 0/(r) = 5,-o and exponential functions 0.s(r) =e^°'^'^, and 27 
choices of p^'^ with Ip^'^l <V3, which cover a region of the momentum transfer — 2.0<^^[GeV^] < 
0. Note that meson momenta p^'' are shown in units of 271/ (A^sfl) in this article. 



3. Electromagnetic form factors 

We calculate an effective value of the pion EM form factor from 

''^^'^^^'^^''^-q;(A.;p)q;,(A.';pO' ^'-'^ 



where we use the dispersion relation to estimate £';i:(|p|) = + p^ at |p| 7^0. As shown in 
the left panel of Fig. |T] we can reliably identify the plateau of Fy^{At,At';q^) by using different 
combination of the smearing functions (0,0'). The EM form factor Fy^{q^) is determined by a 
constant fit to Fy^ {At ,At' ;q^). The statistical accuracy is typically a few % because of the average 
over the source location mentioned above. Although we ignore the finite volume correction to 
Fy^{q^) in this preliminary report, it turned out to be comparable with the statistical error of a few 
% level in our previous study in two-flavor QCD at similar values of m„j and NgQ |^. 

As seen in the right panel of Fig. the q^ dependence of Fy^{q^) is close to the p meson pole 
1/(1— ^^/Mp) expected from the vector meson dominance (VMD) hypothesis. We then assume 
that the small deviation due to the higher poles or cuts can be approximated by a polynomial of q^, 
and use a parametrization 

F^{q^) = -r:^7]j^+^l'?' + ^2(<?Y = l + ^(r'>f -?' + •••• (3.3) 
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to determine the pion charge radius {r^)y^ ■ Since the deviation of Fy^ (q^) from VMD is small, we 
obtain reasonable ;^^/d.o.f. ~ 1, and {r^)y^ does not change significantly by the inclusion of higher 
order corrections of 0{q^). 



At NLO of SU{?>) ChPT, (r^)^ is given by 



-3 + 24NL'g)-2Vj,-VK, Vp 



In 



Mi 



{P = K,K), 



(3.4) 
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where A'^= (47r)^ and the renormalization scale 
/I is set to 4kFo. We fix Fq to 52 MeV de- 
termined from our study of the meson decay 
constants This is significantly smaller 
than the phenomenological value ~ 88 MeV 
[^] and enhances the chiral logarithms Vji^k- 
As shown in Fig. ^ the NLO formula fails 
to reproduce the quark mass dependence of 
our data (dashed line). The extrapolation be- 
comes consistent with experiment by includ- 
ing a NNLO analytic term oc m| (solid line). 
We note that significant NNLO contributions 
have been observed also in our two-flavor stud- 
ies in a similar region of /m„j |^. 

The kaon EM form factors are calculated from the ratios (0) and (3^) but with kaon source 
and sink. Results are plotted as a function of in Fig. ||. The neutral kaon form factor F^{q^) 
originates from a difference between the contributions from the down (dj^d) and strange quark 
currents (syns), and is much smaller than the charged one Fy^{q^). We obtain significant signal 
for F^\q^) with our statistical accuracy improved by using the all-to-all propagator. Similar to 
F^^{q^), the q^ dependence of both F^^{q^) and Ff{q^) is close to that of VMD 



Figure 2: Chiral extrapolation of (r ) " using the one- 



loop formula (3.4) (dashed line) and that with a higher 
order analytic term (solid line). The experimental value 
[R] is plotted by a star. 



1 



3 1-^2/^2 3 1-^VM2' 



1 



3 1-^VM2 3 1-^Vm2- 



(3.5) 
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Figure 4: Chiral extrapolation of charge radii of charged (left panel) and neutral kaons (right panel). Crosses 
and stars are experimental values [0] . 



We determine charge radii, (r^)f^ and {r^)y , using a fitting form with these vector meson poles 
plus a polynomial correction up to 0{q^). 

Figure § shows the chiral extrapolation of (r^)y^ and (r^)^ based on NLO ChPT 

^'^)v = 7;^{-^ + ^^NL',)-v^-2vk, {?)^ = v^-Vk. (3.6) 
2NFq 

For , we again observe that the NLO fit leads to a large value of ;^^^/d.o.f. ~ 3.2. The extrapola- 
tion becomes closer to experiment with acceptable ;^^/d.o.f. (~ 0.6) by including a NNLO analytic 
term. 

Since does not directly couple to photons, the NLO expression of (r^)f° does not have 
analytic terms and hence 0{p'^) LECs. The dashed line in the right panel of Fig. ^ is a parameter- 
free prediction with Fq determined from the decay constants. Our data are consistent with this 
NLO prediction. For a more rigorous comparison, calculations with twisted boundary conditions 
are currently underway to reduce the large systematic uncertainty due to the lack of data near q^=0. 



4. Pion scalar form factor 

In this report, we consider the scalar form factor normalized at a reference value of q^^p namely 
^5^(^^)/^5^(^ref)' since it has sufficient information to determine the scalar radius and does 
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not need a non-perturbative renormalization of S. We set l^^^^l to our smallest non-zero \q'^\, where 
we do not need to subtract the contribution from the vacuum expectation value of S. In the left 
panel of Fig. ^, we plot [q^ ) / (^jef ) determined from ratios similar to Eqs. ( 3J_ ) and ( 3^ ) (see 
Ref. 1^ for details). Due to the lack of the knowledge of scalar resonances at the simulated quark 
masses, we parametrize the dependence of F^[q^) by a generic polynomial form 



Fs{q') = FsiO) (^1 + -{rYs q'+d2{q'Y+d,{qy+d,{q 
We then fit results for to the NLO chiral expansion 



(4.1) 



NFi 



{-S + 24N{2LI + L'^)}-Uvn-3VK. 



(4.2) 



The pion-loop logarithm is 6 times larger than that in and is further enhanced by our small 

value of Fq. As shown in the right panel of Fig.||, the NLO expression has a strong m„j dependence 
and can not describe our data leading to ;^^/d.o.f. ~ 0(100). This is largely reduced to ~ 7 by 
including a NNLO analytic term oc suggesting that the consistency with SU (3) ChPT should 
be studied by including full NNLO corrections as in our previous study in two-flavor QCD [0]. 
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5. Kaon weak decay form factors 

We calculate the vector and scalar form 
factors of the K^n decays, namely f+{q^) 
and fo{q^)=f^{q^) + f^q^) q^ [M], - , 
by using ratios of kaon and pion correlators 
proposed in previous studies [^. For instance, 
(<7^) =/- {q^)/ f+ (<?^), which is needed to con- 
vert f+{q^) to fo{q^) (and vice versa), is deter- 
mined from a double ratio 

c^r(A?,A?',p,p')c^;f^(A?,A?',p,p') ■ 

This involves the three-point functions with spatial components Vk and non-zero meson momenta 
p^'\ that are quite noisy if naively calculated. By using the all-to-all propagator, we obtain a clear 
signal for ^{q'^) as shown in Fig. ^ Both of f+{q^) and fo{q^) are determined with the statistical 
accuracy of typically a few % even at nonzero p(''. 

We plot /() [q^) as a function of q^ in the left panel of Fig.^ The q^ dependence of both /+ o (q^) 
is well described by either the single pole or quadratic form 



Figure 6: Effective value of ^ extracted from ratio 



Eq.(5.1) 



fx 



fx{q') = l+cxq' + dx{q^f (X = +,0) 



l-q^/MY 



(5.2) 



The normalized slopes —^^^hy^^x are measured in recent experiments |^]. Our results = 
2.01(25) X 10^^ and Aq = 1.54(20) x 10^^ extrapolated to M^r phys are in good agreement with the 
experiments as shown in the right panel of Fig. ^. We also observe that the normalized curvature 
A" = 2M^p^y^(i+ = 0.08(10) X 10"^ is also consistent with the experimental value 0.20(5) x 10^^ 
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Figure 7: Left panel: fo{q^) as a function of q^. The solid line shows the quadratic parametrization in 
Eq. (x2). Right panel: linear chiral extrapolation of Aq compared with experiment [Q]. 



6. Conclusions 

We report on our calculation of the Light meson form factors in three-flavor QCD with overlap 
quarks. For the EM and scalar form factors, we observe that the mild rriud dependence of our data 
can not be described by NLO ChPT. We are planning to extend our analysis to NNLO. To this end, 
it is helpful to calculate various observables in order to constrain many and 0{p^) LECs 

involved in NNLO chiral expansions. For instance, we calculate the kaon EM form factors with 
small additional costs using the all-to-all propagator. 

We also confirm that the shape of the K^n decay form factor is in good agreement with ex- 
periment. Our calculations are being extended to a larger volume with twisted boundary conditions 
in order to carry out a chiral extrapolation of /+ (0) with controlled systematic uncertainties, which 
is essential for a reliable estimate of \Vus\. 

Numerical simulations are performed on Hitachi SRI 1000 and IBM System Blue Gene Solu- 
tion at High Energy Accelerator Research Organization (KEK) under a support of its Large Scale 
Simulation Program (No. 09/10-09). This work is supported in part by the Grant-in-Aid of the 
Ministry of Education, Culture, Sports, Science and Technology (No. 21674002, 21684013 and 
220340047) and the Grant-in-Aid for Scientific Research on Innovative Areas (No. 20105001, 
20105002, 20105003 and 20105005) 
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